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ABSTRACT: In this paper we introduce an approach for prognosis of mass transport in heterostructures during 
manufacturing of circuits of a downconversion mixer. The approach gives a possibility to take into account spatial 
and temporal variations of parameters of mass transport at one time. At the same time the approach gives a 
possibility to take into account nonlinearity of mass transport and analyze this process without cross linking of 
solutions on interfaces of heterostructures. Based on this approach, we formulate recommendations for 
optimization of technological process to decrease dimensions of elements of considered circuits and to increase 
their density. Using the approach gives us possibility to take into account mismatch-induced stress, which presents 
in any heterostructures, and porosity of materials of heterostructure. The porosity in our case gives a possibility to 


decrease the considered stress. 
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1. INTRODUCTION 

In the present time several actual problems of the 
solid state electronics (such as increasing of 
performance, reliability and density of elements of 
integrated circuits: diodes, field-effect and bipolar 
transistors) are intensively solving [1-6]. To increase 
the performance of these devices it is attracted an 
interest determination of materials with higher values 
of charge carriers mobility [7-10]. One way to decrease 
dimensions of elements of integrated circuits is 
manufacturing them in thin film heterostructures [3- 
5,11]. In this case it is possible to use inhomogeneity of 
heterostructure and necessary optimization of doping 
of electronic materials [12] and development of 
epitaxial technology to improve these materials 
(including analysis of mismatch induced stress) [13- 
15]. An alternative approaches is to increase the 
dimensions of integrated circuits by using of laser and 
microwave types of annealing [16-18]. 


In the present paper we introduce an approach to 
manufacture field-effect transistors. The approach 
gives a possibility to decrease their dimensions with 
increasing their density framework a down-conversion 
mixer circuit. We also consider possibility to decrease 
mismatch-induced stress to decrease quantity of 
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defects, generated due to the stress. Also, we consider a 
heterostructure, which consist of a substrate and an 
epitaxial layer (see Fig. 1). We also consider a buffer 
layer between the substrate and the epitaxial layer. The 
epitaxial layer includes into itself several sections, 
which were manufactured by using other materials. 
These sections have been doped by diffusion or ion 
implantation to manufacture the required types of 
conductivity (p or n). These areas became sources, 
drains and gates (see Fig. 1). After this doping it is 
required annealing of dopant and/or radiation defects. 
The aim of this present work is to analysis of 
redistribution of dopant and radiation defects and to 
determine conditions, which correspond to decreasing 
of elements of the considered voltage reference and at 
the same time to increase their density. At the same 
time we consider a possibility to decrease mismatch- 
induced stress by diffusion or ion implantation to 
manufacture the required types of conductivity (p or 
n). These areas became sources, drains and gates (see 
Fig. 1). After this doping it is required annealing of 
dopant and/or radiation defects. The aim of this 
present work is to analysis of redistribution of dopant 
and radiation defects and to determine conditions, 
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which correspond to decreasing of elements of the 
considered voltage reference and at the same time to 
increase their density. At the same time we consider a 
possibility to decrease mismatch-induced stress 


Method of solution 
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To solve our aim we determine and analyzed spatio- 
temporal distribution of concentration of dopant in the 
considered heterostructure. We determine the 
distribution by solving the second Fick's law in the 
following form [1,20-23] 


Fig. 1a. Structure of the considered mixer [19] 


Epitaxial layer 


Butter layer 


Substrate 


Fig. 1b. Heterostructure with a substrate, epitaxial layers and buffer layer (view from side) 


AC(x,y,z,t) 0 


polls y,Z,t) i O 


0 


pocluy. zt), 4 


Ot Ox Ox Oy 


oO po luy.zt) 
Oy Ox Oz 


+ Q— De ia (e mee ieee y,W,t)dW |+ 
0 


Ox| kT 


L, 
so! Ps us (x,y,z,t)fC(x,y,W.t)dW + (1) 
Oy| kT 0 
Pankratov 


International Journal of Advanced Science and Engineering 


www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.7 No.2 1723-1751 (2020) 1725 


E-ISSN: 2349 5359; P-ISSN: 2454-9967 


Des O L(x, y, Zt) 


O Des O u(x, y,z,t) es O 
Ox|VkT Ox 


with boundary and initial conditions 


aC(x, y,Z,t) _0 éC(x, y, Zt) 


Ox Ox 


x=0 x=L, 


aC (x, y,z,t) 
Oy OZ 


x=L, 
Here C(x,y,z,t) is the spatio-temporal distribution of 
concentration of dopant; Q is the atomic volume of 
dopant; Vs is the symbol of surficial gradient; 


iL, 
[C(x, y, z,t)d Z is the surficial concentration of 
0 


dopant on interface between layers of heterostructure 
(in this situation we assume, that Z-axis is 
perpendicular to interface between layers of 
heterostructure); su(x,,z,t) and so(x,y,z,t) are the 


C’ (x, y,z,f) 

D.=D,(x, y,z,T ) 1 > 
PMT Br oT) 
D.=D (x, y,z,T) 1+é Ca): 
S SL S P’ (x, y, 27) 


Here Dz (x,y,z,T) and Ds (x,y,z,T) are the spatial (due to 
accounting all layers of heterostruicture) and 
temperature (due to Arrhenius law) dependences of 
dopant diffusion coefficients; T is the temperature of 
annealing; P(x,y,z,T) is the limit of solubility of dopant; 
parameter vy depends on properties of materials and 
could be integer in the following interval ye[1,3] [24]; 


V (xy,z,t) is the spatio-temporal distribution of 
OI(x, y,Z,t) O OI (x, y, Zt) 
or D 9 Noa“ T — 
Bt Ba Hat 
Ol Go y,Z, t) 


O 
oa D a) al 
ae 1 9.2 ) Oz 


Des O u(x, y,z,t) +4 O 
Oy|VkT 


= 1): 


AC (x, ie z,t) 


Oy Oz|VkT OZ 


éC(x, y, Zt) 


=0,c (xy,z,0)=fc (x,y,z), 
Oy 


y=0 
AC (x, y, z,t) 
Oz 


=0, =0. 


z=0 x=L, 


chemical potential due to the presence of mismatch- 
induced stress and porosity of material; D and Ds are 
the coefficients of volumetric and surficial diffusions. 
The values of dopant diffusions coefficients depends on 
properties of materials of heterostructure, speed of 
heating and cooling of materials during annealing and 
spatio-temporal distribution of concentration of 
dopant. Dependences of dopant diffusions coefficients 
on parameters could be approximated by the following 
relations [24-26] 
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concentration of radiation vacancies; V“ is the 
equilibrium distribution of vacancies. Concentrational 
dependence of dopant diffusion coefficient has been 
described in details in [24]. Spatio-temporal 
distributions of concentration of point radiation 
defects have been determined by solving the following 
system of equations [20-23,25,26] 
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Here / (x,y,z,t) is the spatio-temporal distribution of | equilibrium distribution of interstitials; 

concentration of radiation interstitials; /* is the Dv(x,y,Z,T), Dis(x,y, Z,T), Dvs(x,y,z,T) are the coefficients 
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of volumetric and surficial diffusions of interstitials and 
vacancies, respectively; terms V2(x,y,z,t) and /?(x,y,z,t) 
correspond to generation of divacancies and 
diinterstitials, respectively (see, for example, [26] and 
appropriate references in this book); kziv(x,y,z,T), 
kii(xy,z,T) and kvyv(xy,z,T) are the parameters of 
recombination of point radiation defects and 
generation of their complexes; k is the Boltzmann 


constant; w = a3, a is the interatomic distance; l isthe 
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layers we assume, that porous are approximately 


cylindrical with average values f = ae =p y; and 


zi before annealing [23]. Smaller pores is transforming 
to single vacancies during annealing. Single vacancies 
could be absorbed by larger pores [27]. With time large 
pores became larger due to absorbing the vacancies 
and became more spherical [27]. Distribution of 
concentration of vacancies in heterostructure, existing 


due to porosity, could be determined by summing on 


all pores, i.e. 
=x ty tz. 
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and diinterstitials & (x,y,z, t) could be determined by 
solving the following system of equations [25,26] 
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Here Dai(x,y,z,T), Dov(x,y,Z,T), Dais (X,y,z,T) and Devs(x,y,z,T) are the coefficients of volumetric and surficial diffusions 
of complexes of radiation defects; ki(x,y,z,T) and kv(x,y,z,T) are the parameters of decay of complexes of radiation 
defects. 


Chemical potential su in Eq.(1) could be determine by the following relation [20] 
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Farther we replace concentrations of dopant and radiation defects in right sides of Eqs. (1a), (3a) and (5a) on their 
not yet known the average values ai, In this situation we obtain equations for the first-order approximations of the 
required concentrations in the following form 
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We determine average values of the first-order approximations of concentrations of dopant and radiation defects 
by the following standard relation [28] 
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Average values of the second-order approximations of required approximations by using the following standard 
relation [28] 
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Substitution of the relations (1e), (3e), (5e) into relation (10) gives us possibility to obtain relations for required 
average values a2p 
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first-order approximations of the considered known average values ai. The substitution leads to the 
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Approximations of the second and higher orders of components of displacement vector could be determined by 
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Integration of the left and right sides of the above relations on time t leads to the following result 
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Framework this paper we determine concentration qualitative analysis and to obtain some quantitative 
of dopant, concentrations of radiation defects and results. All obtained results have been checked by 
components of displacement vector by using the comparison with results of numerical simulations. 
second-order approximation framework method of 
averaging of function corrections. This approximation 
is usually enough good approximation to make 
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DISSCUSSION 


In this section we analyzed dynamics _ of 
redistributions of dopant and radiation defects during 
annealing and under influence of mismatch-induced 
stress and modification of porosity. Typical 
distributions of concentrations of dopant in 
heterostructures are presented on Figs. 2 and 3 for 
diffusion and ion types of doping, respectively. These 
distributions have been calculated for the case, when 
value of dopant diffusion coefficient in doped area is 
larger, than in nearest areas. The figures show, that 
inhomogeneity of heterostructure gives us possibility 
to increase compactness of concentrations of dopants 


C(x.0) 


0 L/4 
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and at the same time to increase homogeneity of 
dopant distribution in doped part of epitaxial layer. 
However framework this approach of manufacturing of 
bipolar transistor it is necessary to optimize annealing 
of dopant and/or radiation defects. Reason of this 
optimization is following. If annealing time is small, the 
dopant did not achieve any interfaces between 
materials of heterostructure. In this situation one 
cannot find any modifications of distribution of 
concentration of dopant. If annealing time is large, 
distribution of concentration of dopant is too 
homogenous. We optimize annealing time framework 
recently introduces approach [29-37]. 


Substrate 


L/2 3L/4 L 
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Fig.2. Distributions of concentration of infused dopant in heterostructure from Fig. 1 in direction, which is 
perpendicular to interface between epitaxial layer substrate. Increasing of number of curve corresponds to 
increasing of difference between values of dopant diffusion coefficient in layers of heterostructure under condition, 
when value of dopant diffusion coefficient in epitaxial layer is larger, than value of dopant diffusion coefficient in 


substrate 
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Fig.3. Distributions of concentration of implanted dopant in heterostructure from Fig. 1 in direction, which is 
perpendicular to interface between epitaxial layer substrate. Curves 1 and 3 corresponds to annealing time © = 
0.0048 (L,?+Ly*+L2*)/Do. Curves 2 and 4 corresponds to annealing time © = 0.0057(L,?+Ly*+L2*)/Do. Curves 1 and 2 
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correspond to homogenous sample. Curves 3 and 4 corresponds to heterostructure under condition, when value of 
dopant diffusion coefficient in epitaxial layer is larger, than value of dopant diffusion coefficient in substrate 


optimal values of annealing time by minimization of 


Framework this criterion we approximate real 
the following mean-squared error 


distribution of concentration of dopant by step-wise 
function (see Figs. 4 and a Farther we determine 

L,LyL, 
—— || j[C(x.y,z,0)-y(xy,<)] ded ydx, (15) 
i LL, 000 
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Fig. 4. Spatial distributions of dopant in heterostructure after dopant infusion. Curve 1 is idealized distribution of 
dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. Increasing of number of 


curve corresponds to increasing of annealing time 
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Fig. 5. Spatial distributions of dopant in heterostructure after ion implantation. Curve 1 is idealized distribution of 
dopant. Curves 2-4 are real distributions of dopant for different values of annealing time. Increasing of number of 


curve corresponds to increasing of annealing time 
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where wy (x,y,z) is the approximation function. 
Dependences of optimal values of annealing time on 
parameters are presented on Figs. 6 and 7 for diffusion 
and ion types of doping, respectively. It should be 
noted, that it is necessary to anneal radiation defects 
after ion implantation. One could find spreading of 
concentration of distribution of dopant during this 
annealing. In the ideal case distribution of dopant 
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achieves appropriate interfaces between materials of 
heterostructure during annealing of radiation defects. 
If dopant did not achieve any interfaces during 
annealing of radiation defects, it is practicably to 
additionally anneal the dopant. In this situation 
optimal value of additional annealing time of implanted 
dopant is smaller, than annealing time of infused 
dopant. 


@D,L 
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Fig.6. Dependences of dimensionless optimal annealing time for doping by diffusion, which have been obtained by 
minimization of mean-squared error, on several parameters. Curve 1 is the dependence of dimensionless optimal 
annealing time on the relation a/L and é= y= 0 for equal to each other values of dopant diffusion coefficient in all 
parts of heterostructure. Curve 2 is the dependence of dimensionless optimal annealing time on value of parameter 
€ for a/L=1/2 and € = y= 0. Curve 3 is the dependence of dimensionless optimal annealing time on value of 
parameter é for a/L=1/2 and ¢= y= 0. Curve 4 is the dependence of dimensionless optimal annealing time on value 
of parameter yfor a/L=1/2 and e= €=0 
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Fig.7. Dependences of dimensionless optimal annealing time for doping by ion implantation, which have been 
obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and € = y= 0 for equal to each other values of dopant 
diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless optimal annealing 
time on value of parameter ¢ for a/L=1/2 and € = y= 0. Curve 3 is the dependence of dimensionless optimal 
annealing time on value of parameter & for a/L=1/2 and ¢ = y= 0. Curve 4 is the dependence of dimensionless 
optimal annealing time on value of parameter yfor a/L=1/2 and e= €=0 
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Zz 


Fig. 8. Normalized dependences of component uz of displacement vector on coordinate z for nonporous (curve 1) 


and porous (curve 2) epitaxial layers 
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Fig. 9. Normalized dependences of vacancy concentrations on coordinate z in unstressed (curve 1) and stressed 


(curve 2) epitaxial layers 


Further, we analyzed influence of relaxation of 
mechanical stress on distribution of dopant in doped 
areas of heterostructure. Under following condition a< 
0 one can find compression of distribution of 
concentration of dopant near interface between 
materials of heterostructure. Contrary (at a>0) one 
can find spreading of distribution of concentration of 
dopant in this area. This changing of distribution of 
concentration of dopant could be at least partially 
compensated by using laser annealing [37]. This type 
of annealing gives us possibility to accelerate diffusion 
of dopant and another process in annealed area due to 
inhomogenous distribution of temperature and 
Arrhenius law. Accounting relaxation of mismatch- 
induced stress in heterostructure could leads to 
changing of optimal values of annealing time. At the 
same time modification of porosity gives us possibility 
to decrease value of mechanical stress. On the one hand 
mismatch-induced stress could be used to increase 
density of elements of integrated circuits. On the other 
hand could leads to generation dislocations of the 
discrepancy. Figs. 8 and 9 show distributions of 
concentration of vacancies in porous materials and 
component of displacement vector, which is 


perpendicular to interface between 


heterostructure. 


CONCLUSION 


In this paper, we proposed model redistribution of 
infused and implanted dopants with account relaxation 
mismatch-induced stress during manufacturing field- 
effect heterotransistors framework a downconversion 
mixer circuit. We formulate recommendations for 
optimization of annealing to decrease dimensions of 
transistors and to increase their density. We formulate 
recommendations to decrease mismatch-induced 
stress. Analytical approach to model diffusion and ion 
types of doping with account concurrent changing of 
parameters in space and time has been introduced. At 
the same time the approach gives us possibility to take 
into account nonlinearity of considered processes. 
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